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Abstract. We explore the secular dynamical evolution of an N-body model 
of M31 in the presence of a population of 100 dark matter satellites over 10 Gyr. 
The satellite population has structural and kinematic characteristics modelled to 
follow the predictions of ACDM cosmological simulations. Vertical disk heating 
is a small effect despite many interactions with the satellite population with only 
a 20% increase in vertical velocity dispersion <j z and the disk scale height Zd at 
the equivalent solar radius R = 2.5Rd- However, the stellar disk is noticeably 
flared after 10 Gyr with Zd nearly doubling at the disk edge. Azimuthal disk 
heating is much larger with <jr and a z both increasing by 1.7x. However, in 
a control experiment without satellites dispersion increases by 1.5 x suggesting 
that most of the effect is due to heating through scattering off of spiral structure 
excited by swing-amplified noise. Surprisingly, direct impacts of satellites on 
the disk can excite spiral structure with a significant amplitude and in some 
cases impacts close to the disk center also induce the bar instability. The large 
number of dark matter satellite impacts expected over a galaxy's lifetime may be 
a significant source of external perturbations for driving disk secular evolution. 



Introduction 



Cosmological simulations in ACDM show that the dark matter h alos of galax- 



ies contain hundreds to thousands of subhalos or dark satellites ( Moore et al 



19991 : iGao et ahl |2004| ; iDiemand et al.l 120071 ) . The existence of subhalos raises 
interesting questions about the dynamical evolution of disk gal axies. While 



this population is under-represented in the observed satellites (jKlvpin et al 



1999), the recent success of ACDM in accounting for many properties of the 
universe from the CMB to the large-scale structure leads us to take this pre- 
diction seriously. Cosmological infall of a large number of dark satellites onto 
a typical spiral galaxy could be a major source of disk heati ng and thickening 
dToth fc Ostrikerlll992l : iFont et al.l [200ll : iBenson et al] 12004 : iKazantzidis et al.1 
20071 ) beyond known astrophysical processes so it important to quantify the effect 
and see if the observed galaxy population is morphologically consistent. A main 
difficulty with numerical studies is that N-body disks are prone to self-heating 



department of Astronomy and Astrophysics, University of Toronto, Toronto, ON M4S 3H4, 
Canada dubinski@astro.utoronto.ca; nickerson.sarah@gmail.com 

2 Department of Astronomy and Astrophysics, Kavli Institute for Cosmological Physics, The 
University of Chicago, Chicago, IL 60637, USA gauthier@oddjob.uchicago.edu 

3 Department of Physics, Queen's University, Kingston, ON K7L 3N6, Canada 
widrow® astro .queensu.ca 

1 



2 



DUBINSKI et al. 



by two-body relaxation with inadequate resolution so that any perturbations by 
satellites can easily be masked by numerical effects. Ad hoc assumptions about 
the structure of satellites also make the interpretation of results difficult in light 
of the predictions of ACDM. 

In this study, we attempt to overcome these problems using new galaxy 
simulations with sufficient resolution to measure the heating directly due to a 
cosmologically inspired model population of dark matter satellites. These simu- 
lations reveal the importance of satellite impacts - i.e direct passage of satellites 
through the stellar disk - in exciting both spiral and bar instabilities. Dark 
satellite interactions may therefor e be an essential d r iver o f secular dynamical 
evolution of disk galaxies (see also Kazantzidis et all ( 20071 )). 



Methods 



For the study here, we use the stable, equilibrium axisymmetric model of M3 1 
derived from a distribution function as described in lWidrow Dubins ki (2005). 



We model M31 with an exponential disk, a Hernquist model bulge, and a cuspy 
dark halo with an NFW profile. Model parameters are determined that fit the 
rotation curve and surface brightness profile for M31 with an assumed M/L 
ratio for the stellar components such that the disk remains stable against bar 
formation for 10 Gyr. 

The method for generating a satellite population is described in detail in 
IGauthier et al.l (2006). The satellite properties reflect cosmological numerical 



predictions for the subhalo mass function, radial distribution, tidal radii as well 
as internal density structure. To summarize, 10% of the dark matter halo mass is 
initially in 100 subhalos spanning a mass range of 10 8 — 10 10 M selected from a 
mass function dN/dM ~ M . The radial number de nsity of satellites are set 
according to the formulae presented in Gao et all ( 20041 ). Our highest resolution 



galaxy model contains the following numbers of particles: 10M disk, 5M bulge, 
and 20M smooth halo for a total of 35M. The 100 dark satellites each contain 
100K par ticles for a tota l of 10M. We run the simulation using a parallelized 
treecode ( Dubinski 19961 ) with a fixed Plummer softening length e = 15 pc and 



20000 equal timesteps with St = 0.49 Myr. Total binding energy is conserved to 
within 0.3% and total angular momentum is conserved to within 2%. We also 
use models with 10 times fewer particles in some additional studies. 



3. Disk Heating 

We first ran a control simulation at 35M particles to understand numerical ef- 
fects (Fig. 1). Vertical disk heating is negligible over most of the disk with 
almost no change in u z and the vertical scale height at this resolution. However, 
the radial and azimuthal velocity dispersions grow by 50% due to spiral instabil- 
ities arising from swing amplified Poisson noise in the disk particle distribution. 
We then added satellites in two different runs with statistically similar distribu- 
tions. In our first sim ulation, the disk developed a bar half way through the run 
( Gauthier et al.|[20 06) and so introduces an unwanted additional source of heat- 
ing (see Movie 1). The formation of a bar was unexpected and we will discuss the 
origin of this instability shortly. Fortunately in the second run, no bar formed 




Figure 1. The evolution of the disk velocity ellipsoid for the control simu- 
lation (no subhalos) and the simulation with 100 dark satellites with no bar 
instability. Also, the evolution of the vertical scale-height as measured by the 
variance in z at different radii. 



so we could directly measure heating effects by the satellite population (Fig. 1). 
Vertical heating is still small with roughly a 20% increase in a z even under the 
bombardment of satellites. The increase in a T and gq is also only 20% over and 
above the heating caused by intrinsic spiral instabilities. However, there is a 
noticeable flaring of the disk in the presence of satellites wit h the scale height 



nearly doubling from 2 disk scale lengths to the disk edge. iKazantzidis et al 



rly 
)7) 



(2007) also see this effect in similar work. Similar features are seen in the outer 



disk of M31 (llbata et al. 1 12007V ) 



4. Spiral and Bar Instabilities from Satellite Impacts 

An unexpected feature of these simulations is the induced bar instability in 
our first run as well as easily distinguished multi-armed global spiral structure. 
The main cause of these features are satellite impacts. The passage of a satel- 
lite through the disk induces a localized di sturbance that presumably grows by 
Toomre's swing amplification mechanism ( Toomrel 198ll ). Note that the tidal 



effects of the satellites are generally small and so this mechanism is quite dif- 
ferent than the tidal interactions responsible for grand-design spiral galaxies 
like M51. The mechanism is close r to the original i deas suggested by both 
Goldreich Sz Lvnden-Bell ( 19651 ) and I Julian Toomre ( 19661 ) where a mass per- 



turbation appearing within the disk - a giant molecular cloud or massive star 
forming region - is the source of a disturbance that is subsequently amplified. 

A virtual fly by of the evolving galaxy clearly shows episodes of spiral struc- 
ture excitation immediately after satellites pass through the disk (see Movie 2). 
In the case of model with the bar instability, there appears to be a single strong 
encounter with one of the more massive satellites near the center of disk just 
before the onset of the bar. It seems likely the disturbance caused by the passing 
satellite disrupts the center of the galaxy en ough to make it susceptible to the 
swing- amplifier feedback loop ( Toomrelll98ll ). 



We are testing these ideas further with more experiments and idealized 
perturbations representing satellite passages throught the disk. We performed 10 
additional experiments with 4.5M particle models and statistically independent 
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but consistent satellite populations in orbit around M31 (see Movie 3). Five out 
of ten of these models develop a bar instability during their lifetime apparently 
due to chance central encounter with a massive satellite. Those models that do 
not suffer strong central encounters with satellites avoid bar formation. In all 
cases, spiral perturbations of significant amplitude are observed due to satellite 
interactions. 

We are now doing a quantitative study of the effect of satellite impacts 
on disks using a transient mass perturbation appearing within the plane of the 
disk. Preliminary results suggest that even satellite masses as small as 10 9 M Q 
or roughly 5-10% the mass of the LMC are large enough to induce obvious spiral 
structure while the passage of a LMC-sized satellite creates a strong response 
(see Movie 4). A typical galaxy will experience dozens of impacts with satel- 
lites more massive than 10 9 Mq during its life according to the predictions of 
ACDM. Dark satellite impacts may then play an important role in maintaining 
the multi-armed, global spiral patterns seen in the disk galaxies. At this stage, 
we need to quantify the response of a disk as a function of satellite impactor 
mass and impact radius on the disk. The orbital statistics of ACDM subhalos 
from cosmological simulations will allow us to determine the frequency, mass 
distribution and distribution of impact radii expected on a galactic disk over 
its history. By combining these two results, we should be able to quantify the 
effect of dark satellites on disk secular evolution and so address observations of 



the m orphological appearance of galaxies including the bar fraction (jJogee et al 
20041 ) throughout cosmic history. 
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MOVIES 1,2,3 & 4 are available at: 



|http://www.cita.u toronto.ca/ ~dubinski/Rome2007 
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